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Abstract 

We focus on the ground state phase diagram of a system of spinless fermions with repulsion 



> 

on a hypercubic lattice in the limit of infinite dimensions. It occurs spontaneous symmetry 
breaking into a charge density wave (CDW) . Using an ansatz for the order parameter which 
■^j- ' includes the homogeneous, the AB- and a large class of incommensurate phases we are able 

to calculate the phase diagram. 

B 

It is an important object to disclose the structure of the phase diagram of models with strong 
electron correlation. The problem is twofold. First, one has to have a clear picture which phases 
are likely to be realised and second, for these phases it must be possible to calculate the free 
energy. In case of the Hubbard model for example both problems can hardly be surmounted. 
Even if one restricts oneself to a small number of phases and to the dynamic mean field theory 
resulting from the limit of high dimensions [[jj , the problem still requires massive computational 
effort (e.g. pi). The situation is different for the model of interacting spinless fermions H on 
an infinite dimensional hypercubic lattice. Although the diagrammatic theory is simple, the 
physics is far from trivial. 

The Hamiltonian of the spinless fermion model in second quantisation is 

Scaling with the inverse dimension 1/d is performed to ensure the continuity of the limit of 
infinite dimensions |H, 0, . 

In a previous work the stability of the homogeneous phase was investigated in the limit 
d — ► oo [f|. In this limit Hartree- and random phase approximation become exact [||, H. Besides 
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the occurrence of an AB-charge density wave (CDW) characterised by the wave vector Q := 
(7r, 7T, . . . ,7r)t it was found that the density-density susceptibility x(q) = (xo + Ur)(q)j 
diverges at incommensurate values of the wave vector q. The parameter r\ is defined as r](q) := 
1 cos(gi) and Xo l {q) 1S the inverse bare susceptibility. Commensurate values of q which 
belong to a value of r\ ^ —1, 1 are of measure zero. The complete stability analysis is possible 
in d = oo because xo depends on q only via r] Thus it is sufficient to examine xi 7 ]) 0- 

The divergence of x(v) indicates a second order transition from the homogeneous phase at 
high doping and/or low interaction into an incommensurate CDW. The structure of the new 
incommensurate phase has not yet been investigated. In finite dimensions and for related models 
such as the Hubbard model and the t-J model there exist many works which treat different 
variational or approximative approaches to incommensurate phases ||. Since the perturbation 
theory of the model (1) becomes exactly tractable in d — > oo it is the natural candidate to 
examine the energetic effects of different spatial structures of the order parameter. The resulting 
equations are self-consistent and include infinitely many parameters in the thermodynamic limit 
so that no systematic approach to their solution exists || . Here we give results for the following 
ansatz for the spatial structure of the order parameter b(r) := (fir) — n where n is the average 
particle density: b(r) = bo Yli=i u i( r i)- The index i counts the directions; is the iih component 
of r, bo is the overall amplitude and the new defined functions Ui take the values ±1. The product 
form of this ansatz allows to profit best from the simplifications of d — ► oo. Yet it is still very 
general since it allows any sequence of +1 and —1 independently for the d directions without 
assuming translational invariance. 

A close inspection of the free energy shows that it depends on b(r) only via the amplitude 
bo and the relative frequency h that b(r) has the same sign on adjacent sites §|. These are the 
only parameters which are relevant for the influence of the spatial structure on the energy. In 
order to find an explicit functional for the energy we exploit the freedom in choosing the Uj. For 
i < hd we set Ui(ri) = 1 and for i > hd we set Ui(ri) = (— l) r \ This choice is special since the 
order parameter can be characterised by the wave vector 

Q h = (0,0 Q,ir,f ,i)t (2) 

hd {l~h)d 

but does not restrict the freedom contained in the ansatz for b(r). This gives a general relation- 
ship r\ = 2h — 1, although only for special realisations r/ can be defined. The easy structure of 
the CDW with wave vector Q h helps to resum the perturbation series for the grand canonical 
potential since the Green function is a 2 x 2 matrix in fc-space. 

In infinite dimensions the diagram for the free energy is simple. Compared to the free fermion 
system the grand potential gains a term Aft = TTrf^ dXX' 1 GE| u^ux, where G is the Green 
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function and S the self-energy. At T = this yields the ground state energy 
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P(A, 77) = — — - — — + dv — , \ '{ I -— erf (P) - v ' '- exp(-P 2 ) } (3) 
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where \x = [i — nil and A = —r/Ubo. The doping <5 = 0.5 — n is given by 
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dv 9 7 n ^ erf ( p ) ( 4 ) 
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and the equations which define the energetic minimum are dE/dA\ n = and dE/drj\ n = 0. For 
a particular value of the interaction U it is possible to determine for each doping the parameters 
A and rj which describe the physical state. Additionally, one has to be on the watch for first 
order phase transitions. 

As example we look at the case U = 1.2. We plotted ground state energy as function of 
doping in Fig. 1. For large doping and at half-filling the situation is clear. For the former the 
system is in the homogeneous phase whereas for the latter the system is in the commensurate 
CDW characterised by 77 = —1 (AB-phase). Inbetween four special points at 5 div _ AB , <WiP) 
<5ps-ab, and <5 div _ IP can be identified. 

Let us focus on the situation with fixed 77 = — 1. The system is forced to choose between the 
homogeneous phase and the AB-phase. The response of the system towards fluctuations with 
wave vector Q becomes infinite at cW-ab- The ground state energy in indicated by a dotted 
and a dotted-dashed curve, respectively in Fig. 1. For smaller doping the ground state energy 
is concave so that a first order phase transition will occur. The Maxwell construction gives rise 
to the value cWab- Thus, there is phase separation between the homogeneous phase and the 
AB phase for < 5 < 5 ps .ab- This situation turns out to be generic for U < U 1PL ~ 0.5716 and 
U > U IPH ~ 1.9145. Inbetween - and thus for the case U = 1.2 - a different scenario is realised. 

The divergence of x(v) is a t ^div-ip for general values of 77. For ^7 IPL < U < U IPH holds 
<Jdiv-ip > <5ps-ab- By minimizing E(A,rj) we obtain the solid curve in Fig. 1. Again, the curve 
is concave for low dopings so that a Maxwell construction must be made. The doping where 
the phase separation sets in is <5 PS _ IP . It is important to realise that S PS _ IP < <Wab- F° r dopings 
< 5 < <5 div _ IP the pure incommensurate phase is the absolute energy minimum. 

The overall result is depicted in Fig. 2. At half-filling the system is in the AB-phase for all 
U > 0. For large doping the system is in the homogeneous phase. For U < f7 IPL and U > f7 IPH 
there is phase separation between the AB-phase and a hole-enriched homogeneous phase. For 
U lPL < U < Uiph there exists a region where an incommensurate phase is present, and a region 
of phase separation between the AB-phase and the hole-enriched incommensurate phase. 
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The parameter rj is close to —1 in the incommensurate phase. To give an impression we 
plotted those values of r\ in Fig. 3 which arise at the borders of the incommensurate phase in 
Fig. 2. The function r](8) at constant U is monotonic so that the figure gives a good impression 
of the deviation of rj from —1. Note that the values for r/ belonging to the left border (Fig. 2), 
given by the solid curve in Fig. 3, are important for the whole region of phase separation. 

In summary, we demonstrated the complexity of the ground state phase diagram of the 
model of interacting spinless fermions in infinite dimensions. Most important is the existence of 
an incommensurate phase and of the phase separation between the AB-phase and the incom- 
mensurate or homogeneous phase. 

This work was supported by the Deutsche Forschungsgemeinschaft (SFB 341) and the Sticht- 
ing voor Fundamenteel Onderzoek der Materie (FOM). 
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Figure captions 

Figure 1. Ground state energy for U = 1.2. Dotted curve: homogeneous phase. 
Dotted-dashed curve: AB phase. Solid curve: incommensurate phase (77 optimised). 
Divergence of x(~ 1) a t ^div-AB- Divergence of x(v) f° r optimised 77 at <5 div _ IP . Maxwell 
constructions between 5 = and <5 PS _ AB or <5 PS _ IP . Incommensurate phase between 
<5 PS _ IP and <5 div _ IP . To facilitate visibility a linear function of 5 has been added: a ~ 
0.5561, b ~ 0.1946. 

Figure 2. Phase diagram. HOM: homogeous phase, AB: AB phase, IP: incommen- 
surate phase, PS-AB: phase separated region between homogeneous and AB phase, 
PS-IP: phase separated region between incommensurate and AB phase. C/ IPL , t/ IPH : 
lower and upper interaction bound of the IP. 

Figure 3. The value of 77 along the phase borders of the incommensurate phase. 
Solid curve: right border; dashed curve: left border. 
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